oo 

o 

O 
(N 

D 
tin 



Supersymmetry and Nonequilibrium Work Relations 

Kirone Mallick,^ Moshe Moshe,^ and Henri Orland"^ 

'Service de Physique Theorique, Centre d'Etudes de Saclay, 91191 Gif-sur-Yvette Cedex, France 

^Department of Physics, Technion - Israel Institute of Technology, Haifa 32000, Israel 
^Service de Physique Theorique, Centre d'Etudes de Saclay, 91191 Gif-sur-Yvette Cedex, France 

(Dated: February 2, 2008) 

We give a field-theoretic proof of the nonequihbrium work relations for a space-dependent field 
with stochastic dynamics. The path integral representation and its symmetries allow us to derive 
Jarzynski's equality. In addition, we derive a set of exact identities that generalize the fiuctuation- 
dissipation relations to far-from-equilibrium situations. These identities are prone to experimental 
verification. Furthermore, we show that supersymmetry invariance of the Langevin equation, which 
is broken when the external potential is time-dependent, is partially restored by adding to the 
action a term which is precisely Jarzynski's work. Jarzynski's equality can also be deduced from 
this supersymmetry. 
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During the last decade a number of exact relations have 
been derived for non-equilibrium processes. The Jarzyn- 
ski equality is one of these remarkable results: it shows 
that the statistical properties of the work performed on 
a system in contact with a heat reservoir at temperature 
kT — P~^ during a non- equilibrium process are related 
to the free energy difference AF between two equilibrium 
states of that system. This identity was derived originally 
using a Hamiltonian formulation [l|] and was extended to 
systems obeying a Langevin equation [2| or a discrete 
Markov equation [3|, lj|. Jarzynski's result has been ver- 
ified on exactly solvable models [5| and by explicit cal- 
culations in kinetic theory of gases [g, [^ . This equality 
has also been used in various single-molecule pulling ex- 
periments [sl, y, llO| to measure folding free energies and 
has been checked against analytical predictions on meso- 
scopic mechanical devices such as a torsion pendulum 
ll|. These experiments are delicate to carry out be- 
cause the mathematical validity of Jarzynski's theorem 
is insured by rare events that occur with a probability 
that typically decreases exponentially with the system 
size (for a review see e.g. [12|). 

In the present work, we derive nonequilibrium work re- 
lations for a field (f>{x, t) representing a coarse-grained or- 
der parameter of a microscopic system. This field evolves 
according to an effective stochastic equation that de- 
pends on the symmetries and conservation laws of the 
system [l3[. We represent the stochastic evolution as 



a path integral |14l . Il5l . Il6j and use the response field 



formalism 17[ to derive the work relations from the in- 
variance of the path integral under certain changes of 
variables. We obtain correlator identities that generalize 
the fluctuation-dissipation relations arbitrarily far from 
equilibrium. These identities can be checked experimen- 
tally in single molecule experiments. Furthermore, by 
introducing auxiliary Grassmannian fields, we interpret 
this invariance as a manifestation of a hidden supersym- 
metry. This supersymmetry is known to be the funda- 



mental invariance property that embodies the principle 
of microscopic reversibility and leads to the fluctuation 
dissipation-theorem and the Onsager reciprocity relations 
for a system at thermal equilibrium [l8|, ll9|, l20|, L 



21|. Here, 



we show that, far from equilibrium, Jarzynski's theorem, 
is also a consequence of an underlying supersymmetry. 
This supersymmetry in turn allows to generalize the fluc- 
tuation dissipation-theorem to far from equilibrium situ- 
ations. 

We consider a scalar field (/)(x, t) defined on d-dimensional 
space with Model A dynamics that describes a system 
with non-conserved order-parameter [13| {e.g., the Ising 
model with Glauber dynamics) : 



dt 



/\7 J 

T„—+ax,t) = -Taf[c^) + ax,t) (1) 



C(x, t) is a Gaussian white noise of zero mean value and 
correlations {C(x,t)C(x' ,t')) = 2TQkTS{t - t')6'^{x - x'), 
T being the temperature. The dynamics is thus governed 
by the time-dependent potential 

U[4){x,t),t]^T[(t){x,t)]- I (fxh{x,t)(t){x,t) (2) 



The potential energy (or Euclidian action) J-[4>\ is, for 
instance, given by 



T[<t>] 



dM^^o0' + ^|V0p+,.o'^4| (3) 



When the external applied field h{x, t) is constant in 
time, the invariant measure associated with Eq. ^ is 
the equilibrium Gibbs-Boltzmann distribution: 

PeM = ^p-^ With Z[P, h]= JV<P e-^"[^l (4) 

We now consider the case where the applied field varies 
with time according to a well-defined protocol: For < < 0, 



we have h{x,0) — ho{x) and the system is in its station- 
ary state; for i > 0, the external field varies with time, 
reaches its final value hf{x) after a finite time tj, and 
remains constant for t > tf. The values of the poten- 
tial U ioT t < and t > tj are denoted by Wq and Ui, 
respectively. 

The probabihty V {(t)i\(j)o) of observing the field 0i(a;) 
at time tf starting from (j)oix) at time i = is given by 

T'C^i |0o) ^ fvCe-^^ ''"'"''^'S{<t>ix, tf)-^, (x)) 



Writing the second term on the r.h.s. as fdu — SH 



(5) 
The following identity is now substituted in Eq. ([5]) 

1= f V^{x,t)S U{x,t) +To^ ~ C{x,t)\ |detM| 
V^{x,t)V4>{x,t)\detM\ e-/'^'=^''*^W+ro4^-C} (g) 



where 4>{x, t) is the response field and M is the operator 
,, d ^ df{(t>{x,t),t) 



Integrating over the noise ({x,t), one finds [21|.|22| 



{x,0)=rpo{x) 



where the dynamical action S is given by 






To (50 f3' 2 ^02 



(8) 



(9) 



the last term being the Jacobian of M. Wc consider now 
a functional 0[4i] that depends of the values of the field 
(j){x, t) ior < t < tf. The average of ©[0] with respect to 
the stationary initial ensemble and the stochastic evolu- 
tion between times and tf is given by the path integral 

where 

l'(f>{x,tf)=4>i{x) 

J{0,0i,0o}=/ V<i)(x,t)V^{x,t) 

Under a change of the integration variable cj) in Eq. PH)) . 

JU[(l}ix,t),t] 



{x,t) -^ -(j}{x,t)+l3- 



(12) 



S(j){x,t) 
the path integral measure is invariant and E varies as 

E(0, ^, 4>) ^ S(0, -0, 4>) + /3,^^^M^lM (13) 



d'^xdi 



■ 5U[(j){x,t),t] 
5(j}{x, t) 



rfF-oF^ gives 



^Ui[cj>i]-U^[M-Wj[4>llA) 



The derivative dU/dt is related to Jarzynski's work by 



dt^r- = - / d'^xdt h{x,t)(t){x,t) (15) 







at 



the last equality being obtained from Eq. ^ . The change 
of sign of the time derivative (/> in Eq. P^ is compensated 
by the change of variables in the path integral 

((/.(x, t),^{x, t)) -. {^{x, tf - t), ^{x, tf - t)) (16) 

This time-reversal transformation leaves the functional 
measure invariant and restores S to its original form but 
with a time-reversed protocol for the external applied 
field h{x, t) —^ h{x, tf — t). Performing the above change 
of variables p^ and pB]) in Eq. pT|) and using Eqs. HJ 
and 1131), wc find 



X{0, 01, 0o} = e'3("°[^"l-'^i[^il)x{e-'3i^.'a, 0o, Mr (17) 

The subscript R denotes a time-reversed protocol and the 
time-reversed 0[(t>\ is equal to 0[4){x,tf — t)]. Inserting 
this identity in Eqs. (fTOl - [TT|) , gives 

(O) = ^ /p0o(a;)P0i(x)e-^"^[*il X{e-^^^'O,0o,0i}i 






(18) 



where Ai^ is the free energy difference between the final 
and the initial states. Finally, we redefine O as Oe~^'^' . 
Recalling that the work Wj is odd under time-reversal, 
we deduce from the last equation that 

{Oe-^^-')^e-P^^{d)B. (19) 

When O = 1, we obtain Jarzynski's theorem 

(e-'^^O = e-P""" (20) 

Taking O = e'^'^^^^^-', where A is an arbitrary real pa- 
rameter, we derive the following symmetry property 



(e-^^^)=e-'^^^(e(^-'^^'^'')fl, (21) 

and a Laplace transform leads to Crooks relation [sl, lj| : 



Vf{W) 
Vr{-W) 



^ g/3(W-AF) 



(22) 



where Vf and Vr represent the probability distribution 
functions of the work for the forward and the reverse 
processes, respectively. We emphasize that our proof of 
Crooks and Jarzynski identities is based on invariance 
properties of the path integral and does not involve any 



a priori thermodynamic definition of heat and work. The 
expression (llSp for the Jarzynski work appears as a nat- 
ural outcome of this invariance. 

The identity p^ . which is at the core of the work 
fluctuation relations, is valid for any choice of the external 
field protocol. The free energy variation is a function only 
of the extremal values of the applied field at io = and 
t = tf and is independent of the values at intermediate 
times. Functional derivatives of Eq. ([20ll with respect 
to h{x,t) at an intermediate time tg < t < tf, and at 
position X, results in new identities 



{{^ix,t) 






ix,t))e-^'^-') 







(23) 



The 71-th functional derivative of Eq. (PO)) at intermediate 
times ti, . . .tn and positions xi, . . .Xn, gives the identity 



(e- 



-PWj 



Y[{(l){xi,k) 



^^{x,,U)))^Q (24) 



Similarly, the functional derivative of Eq. p9l) leads to 



1 
To 



((^(x,t) - i-0(a;,i))Oe-'5^') = e-^''^ {4>{x,t)d) r 



Eq.dlSD follows by choosing = 0^1 since [231 (0) = 
0. For the special case 0[4)\ = (j){x',t'), we obtain a 
generalization of the fiuctuation-dissipation theorem 



5{(t){x',t')e''f^^^) 



Shi{x,t) 



hi=a 



P{<P{x,t))cb{x',t')e-f''^-') 



-f3AF 



5{^{x\t'))j 



5hi(x, t) 



/n=o 



(26) 



where hi{x,t) is a small perturbation that drives the 
system out of the protocol h{x,t). Note that hi does 
not enter the definition of the Jarzynski work. This new 
fluctuation-dissipation theorem relates out of equilibrium 
response functions to derivatives of correlation functions 
and could be verifled experimentally, for example in sin- 
gle molecule pulling experiments (this corresponds to the 
case where the field (j) does not depend on space). For 
a system at thermodynamic equilibrium with constant 
external field {i.e., Wj — AF = 0) and stationary corre- 
lations, this equation reduces to the standard fiuctuation- 
dissipation relation [20J 

Identities between correlators such as Eqs. (l23t - (|25)) 
suggest the existence of an underlying continuous sym- 
metry of the system. We first extend the integration 
range of the path integral in Eq. pT7)) over the range 
— oo < t < -l-cxD, using the following properties of the 
probability distribution: 



^^-mAM= lini P{^o\M 



T^-c 



1 = jv<P{x,T)P{cj>T\K) 



for T>ti 



(27) 
(28) 



The first property assumes ergodicity and the latter is 
normalization. Inserting into these equations the path 
integral representation, Eq.® , of P(0'|0"), Eq. (fTO|) is 
rewritten as 

{0)= I Vc^{x,t)VUx,t) e-^'^''"^'^^^^'^^^'^0[cf\ (29) 



where 4>{x,t) and 4>{x,t) are integrated with t ranging 
now from — oo to oo. 

To uncover the above mentioned hidden symmetry, in 
addition to the original field (/'(a:,i) and the response 
field 4>{x,t), we introduce two auxiliary anti-commuting 
Grassmannian fields c{x,t) and c{x,t) that allow us to 
express the Jacobian of M, defined in Eq. ([7]), as a func- 
tional integral 19|, l20| . Assuming that O differs from the 
identity only for < i < t/ , the mean value of O in 
Eq. (flO|) can be rewritten as 



(O) 



V(f>V(j)VcVc e~/'i''^'i*^('^'^'=^s) 0j^] 



(30) 



su 



where S(0, </), c, c) = ro0(-f + — - ^) - cMc (31) 

Fo d<p [j 

with M given in Eq. jT]). 

Consider now the infinitesimal transformation that 
mixes ordinary fields with Grassmannian fields: 

5(f){x,t) = c{x,t)e 54>{x,t) = -— c{x,t)e (32) 

To 



5c{x,t) = 5c{x,t) = 4>{x,t) — -—(j){x,t) I e 

V To / 

e being a time-independent infinitesimal Grassmannian 
field. The variation of S((/), 0, c, c) in Ea. ipT]) under the 
transformation of Eq. ([5^ gives 

^i:(0,0,c,c-) = ^-/3|(^)c(.,t).- (33) 



with the total derivative term 






(34) 



If the potential Li is independent of time, the variation 
of S under the supersymmetric transformation ([32|l is 
a total time derivative that does not modify the action. 
The supersymmetry in Ea. (|32p which is in fact a gener- 
alization of supersymmetric quantum mechanics IT], |2l| , 
refiects the time reversal invariance of Model A in ab- 
sence of external field and allows to prove the fiuctuation- 
dissipation theorem [l9|, |20| . 

When the potential l/({(j),t) depends explicitly on time, 
supersymmetry invariance is broken: the last term in 
Eq. (j33[) breaks the invariance. It can be written as 



6^U _ ^s^u ( du 



(35) 



and can be interpreted as the variation of a function. 
Hence, the modified Sj, defined as 



S^ = S + /3|^ (36) 

and obtained by adding the Jarzynski work psp to 
the initial action, is invariant under the supersynimet- 
ric transformation (|32p up to a total derivative term: 
ST, J = ^. The boundary terms at i = ±c>o are, conven- 
tionally, assumed to vanish. 

We now show that the supersymmetric invariance of 
Sj implies the correlator identities (p4|) . Introducing a 
four-component source {H,H,L,L) , we define the gen- 
erating function 



Z{H,H,L,L) 



the potential varies according to a time-dependent proto- 
col. We have shown that the supersymmetry in Eq. ([33l) is 
restored by adding to the action a counter-term which is 
precisely the Jarzynski work /3W[0]j. Furthermore, we 
proved that the associated supersymmetric Ward Iden- 
tity implies Jarzynski's theorem. Supersymmetry en- 
forces the exactness of the quasi-static limit even for pro- 
cesses that have a finite duration and that bring the sys- 
tem arbitrarily far from equilibrium. A hidden supersym- 
metry [25] is also present in classical Hamiltonian systems 
for which Jarzynski's equality was initially proved. We 
finally remark that supersymmetry may also be a useful 
tool when applied to the fluctuation theorem for stochas- 
tic dynamics [2a |. 



exp 



d'^xdt I 



V(j)V(l)VcVc 



,c,c) + H(j> + H(j) + Lc + Lc) 



Making the transformation ([5^ in Z{H, H, L, L) and us- 
ing the supersymmetric invariance of Sj in Eq. (|36p . we 
deduce as in [20| the Ward-Takahashi identity: 

To dtSL \5H TodtSHj SL J 

(37) 
By applying to the Ward-Takahashi identity the operator 

,rf ,N nr 1 I xu/ 4 \ — r:^ xu/ 4. > ) and setting the 

source field H, H, L, L to zero, we obtain Eqs. (|23p and 
(IM)) . This lead to Jarzynski's equality (^0)1 . Replacing 
h{x, t) by h{x, at) for any a > 0, we obtain 



d{e 



-0Wj 



da 

Using Eq.(l23l 

d(e-'^^'')/da 



= / d'^xdt t h{x,at){{(j){x,t) - -^^{x,t))e 



-pw.^ 



-l3Wj 



that can be obtained from ([37|) . we get 
= 0, which means that the value of 



does not depend on a. Hence, this value is the 
same as that of the quasi-static limit a ^ 0, and is given 
by exp(— Ai^). Jarzynski's identity is thus obtained as a 
consequence of supersymmetry. 

The response- field technique in Eqs.([8lfTT|l that we have 
used to derive nonequilibrium work theorems for Model A 
can be extended to multi-component fields and to other 
stochastic models such as model B. It also can be ex- 
tended to systems with correlated noise [2J| replacing 
the Gaussian measure in the RHS of Eq.® by 



VCe 



-i/d'^s dt d'^y dt' C(x,t) A-^{x,t;y,t') Ciy.t') 



(38) 



where A(a;, t; y, t') is the two point correlation function. 
We have obtained correlators identities involving an 
arbitrary field-operator and also a generalization of the 
fluctuation-dissipation relation that remains valid far 
from equilibrium. The supersymmetric invariance of the 
time independent Langevin equation breaks down when 
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